Density-matrix renormalization group study of pairing when electron-electron and 
electron-phonon interactions coexist: Effect of the electronic band structure 
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Density-matrix renormalization group is used to study the pairing when both of electron-electron 
and electron-phonon interactions are strong in the Holstein-Hubbard model at half-filling in a re- 
gion intermediate between the adiabatic (Migdal's) and antiadiabatic limits. We have found: (i) 
the pairing correlation obtained for a one-dimensional system is nearly degenerate with the CDW 
correlation in a region where the phonon-induced attraction is comparable with the electron-electron 
repulsion, but (ii) pairing becomes dominant when we destroy the electron-hole symmetry in a tres- 
tle lattice. This provides an instance in which pairing can arise, in a lattice-structure dependent 
manner, from coexisting electron-electron and electron-phonon interactions. 
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Introduction — The problem of what happens when 
the electron-electron interaction coexists with electron- 
phonon interaction is arousing interests from various 
viewpoints P, Qi S S 111- In fact, there are a num- 
ber of classes of materials where electron-electron (el-el) 
and electron-phonon (el-ph) interactions are both sig- 
nificant. One example is the superconducting doped 
fullerenes0, where the electrons in narrow conduction 
bands are strongly coupled to intra-molecular, high- 
frequency phonons. Theoretically, a fundamental ques- 
tion is: can superconductivity arise in the coexistence 
of the el-el and el-ph interactions? It may seem diffi- 
cult for them to work constructively, since the effective 
electron-electron attraction arising from the el-ph cou- 
pling favors isotropic pairs, while the (spin- fluctuation 
mediated) pairing interaction due to the el-el repulsion 
favors anisotropic pairs. The problem becomes especially 
nontrivial for a region intermediate between the adiabatic 
(Migdal's) and antiadiabatic limits. 

There is a theoretical reason why we have to look into 
the intermediate regime. Superconductivity (SC) has to 
compete with diagonal orders in general, and SC phases 
in fact often arise adjacent to density- wave phases on the 
phase diagram. In the present context, a strong el-ph in- 
teraction will favor a charge density wave (CDW), while 
a strong el-el interaction will favor a spin density wave 
(SDW). A metallic phase has a chance to appear around 
the boundary where CDW gives way to SDW, at which 
the el-el and el-ph interactions are indeed comparable. 
Takada and coworkers have in fact argued that an off- 
site pairing arises in a region between CDW and SDW 
phases, from the enhancement factors in the CDW, SDW 
and SC response functions in a two-site result by assum- 
ing that the two-site system is already close to the infi- 
nite systeml^, and from an exact diagonalization for two- 
and four-site Holstein-Hubbard model at half-filling 8] . 
However, a study for larger systems may be desirable to 
elaborate why and how the pairing occurs when the el- 



el and el-ph interactions are comparable. This region is 
physically interesting as the case of interacting fermions 
and bosons with similar energy scales, but is technically 
challenging as well, since we are not allowed to adopt 
the Migdal's approximation for the phonon energy fko as- 
sumed to be much smaller than the electron energy scale, 
i, nor can we adopt the antiadiabatic limit (Huj ^ t). 

The present Letter focuses on this problem, where 
we adopt the density-matrix renormalization group 
(DMRG)"^ for a one-dimensional Holstein-Hubbard 
model. DMRG can treat interactions non-perturbatively, 
and various types of correlation functions at T = can 
be obtained from the wavefunctions. A crucial interest 
here is how the system is dissimilar (or similar) to purely 
electronic systems, so we scan the three parameters (the 
three axes in the inset of Fig.l) characterizing the system: 
U/t (the electron-electron repulsion U divided hyt), TiLo/t 
(ratio of the phonon and electron energy scales), and 
A = Ig^ jliio (the effective attraction between electrons on 
the same site as defined in the w — > oo limit). As for the 
electronic band filling, we have concentrated here on the 
half-filling as in (l[lfl[llII3> having AgCgolA = K, Rb) 
in mind. 

We shall conclude that, (i) while the correlation func- 
tions obtained here indicate that superconductivity does 
not dominate over, but is nearly degenerate with the 
CDW correlation (a curious similarity with the behav- 
ior in purely electronic systems where the two correla- 
tions should coincide with each other when the system is 
electron-hole symmetric), but (ii) pairing becomes dom- 
inant when we destroy the electron-hole symmetry (in 
a trestle lattice). These occur in the region of interest, 
where the phonon-induced attraction almost cancels the 
electron-electron repulsion, and intermediate between the 
adiabatic and antiadiabatic limits. So the message here is 
that the coexistence of el-el and el-ph interactions can, in 
a manner dependent on the underlying electronic struc- 
ture, give rise to pairing. 
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Method — Inclusion of the phonon degrees of freedom 
makes the DMRG calculation more demanding. How- 
ever, Jeckelmann and White have introduced the pseudo- 
site method J which makes the application of DMRG 
to models with on-site (Einstein) phonons feasible, and 
charge and spin gaps have been calculated for the one- 
dimensional Holstein-Hubbard model [lllll^ . 
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Here creates an electron of spin a on the ith lattice 
site, t is the nearest-neighbor hopping (which we take as 
the unit of energy hereafter), U is the on-site electron- 
electron repulsion, a\ creates a phonon at ith, g is the 
on-site el-ph interaction, and huj is the bare phonon en- 
ergy. In the pseudo-site method, we consider only finite 
numbers of phonons at each site (M = 2^ for 0-boson, 
1-boson, . . . , 2^ — 1-boson states). Sites are expressed in 
terms of fictitious N phonon "pseudo-sites" , on top of an 
electron pseudo-site. The pseudo-sites are taken into the 
system one by one, as in real sites in conventional DMRG. 
Thus the maximum dimension of the Hilbert space con- 
sidered in the DMRG depends only on the maximum 
number of retained states m and not on M . We have 
here retained up to m = 600 states in each block of the 
DMRG calculation for each chain with L = 40 sites with 
M = 2'^ states considered for each site. 

In calculating the correlation functions, we have no- 
ticed a flaw in the conventional infinite-size DMRG algo- 
rithm, in which two new sites are inserted between two 
symmetrical blocks. For the Holstein-Hubbard model, 
the retained states when the electron pseudo-site is added 
tend to significantly deviate from those describing the 
system at later stages. This is because the electrons on 
the added pseudo-sites feel the bare repulsion there, while 
the repulsion on other sites is reduced due to the electron- 
phonon coupling, and this degrades the convergence. So 
we propose [l^ to remedy this by modifying the chemical 
potential for the electrons for the added electron pseudo- 
sites in the calculation for the ground states, so that the 
expectation value of the number of electrons at each of 
the new pseudo-sites equals those in the target state, 
which is unity for the half- filled band assumed here. This 
method, which we call the compensation method, has in- 
deed given lower ground state energies and hence more 
accurate ground states and correlation functions. 

Results — We have calculated the correlation func- 
tion (0|0j) for Oi = nil -I- (charge), (n^i - nii)/2 
(spin), Ci-^Cii (on-site pair), Cif Ci+i| — Ci+^Cii (nearest- 
neighbor spin-singlet), and Ci-\Ci+ii + Ci+i-^Cii (nearest- 
neighbor triplet). Typical result is displayed, on double- 
logarithmic scales, in Fig^for the cases where (A) U > \, 
(B) U = X, and (C) U < \. After the infinite-algorithm 
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FIG. 1: Various correlation functions versus real-space dis- 
tance in the Holstein-Hubbard model at half-filling calculated 
for (A,Lj) = (3.6,5.0) (in units where t = 1, ?i = 1) with 
three values of (7, indicated as A,B,C in a schematic param- 
eter space (top inset). Error bars are smaller than the size 
of each symbol. Only the data for |i — j| =odd are plotted, 
because of an even-odd effect in correlation functions in open- 
boundary condition. Vertical arrows with (a),(b) correspond 
to the scans in Fig. |5Ja)(b). 



DMRG calculation is performed with the compensation 
method, at least three sweeps of the finite-algorithm are 
done. 

We immediately notice the following: if we concen- 
trate on power-law correlations, they are the charge and 
on-site pair correlations for J7 < A, or the spin corre- 
lation for t/ > A, while all of the charge, spin and on- 
site pair correlations decay with power laws for J7 ~ A. 
This is our first key result. It is interesting to compare 
the present correlation functions with an estimation of 
spin and charge gaps in the same system by Fehske et 
al. (1^ They suggested that a charge gap opens and spin 
is gapless when t/ » A, both gapless when [/ ~ A, both 
gapful when U ^ X. The former two agree with the 
present result. As for the latter, we can interpret that 
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here we are looking at a region where U is only moder- 
ately smaller than A and a spin-gapped metal can exist, 
for which CDW and on-site SC have power-law correla- 
tions as far as the Tomonaga-Luttinger physics0| goes. 

So a next key issue is whether the exponents are solely 
determined by J7 — A, which would be the case if the TL 
picture somehow holds everywhere. We have looked at 
its systematic dependence, and Fig. ^ plots typical be- 
haviors. For a smaller A (cx g'^) the dominant correlation 
in Fig.El;a) changes from the CDW/on-site SC to SDW 
almost exactly at the point where U exceeds A. 

To be more precise, it is difficult, for finite systems, 
to distinguish whether the decay of a correlation func- 
tion is power or exponential. So we have checked how 
the result compares with the exponent, calculated in the 
same way, for a purely electronic (i.e., the Hubbard) 
model for the same system size in Fig.|2Id). In the one- 
dimensional Hubbard model, the correlations that de- 
cay with a power law (~ 1/r) are on-site SC and CDW 
for attractive U < (SDW for repulsive U > 0). In 
Fig. Eld), this appears as a crossing of two, continuous 
curves (one of which are doubly-degenerate) at C/ = 0, 
and the result resembles Fig. |2Ia). For finite systems we 
have also to be careful about the efi'ect of discreteness of 
the levels (a gain st spin and charge gaps) on the correla- 
tion function [l^. We have checked that the spin gap for 
(J7,A,w) = (2.0,3.6,5.0) extrapolated from the L < 50 
result is about twice as large as the difference between 
the lowest-unoccupied and highest-occupied one-electron 
levels. So finite-size effects would not be significant, at 
least for \U -X\> 1.0. 

However, we have found that, when A is increased (with 
a fixed adiabaticity Lu/t), the crossing point does deviate 
from [/ — A = as seen in Fig. ^h). Also, the degen- 
eracy between CDW and on-site SC is lifted (with the 
power for the latter slightly larger). This becomes more 
manifest for smaller lu, as shown in Fig.|2Ic) for w = 2.5. 
There again, CDW and on-site SC correlations both have 
power-decay, but they start to decay exponentially when 
U is further decreased, signalling an opening of the charge 
gap Ac as well as spin gap A^. This should correspond 
to the region of A^ - A^ > for [/ < A in Ref.[li|. 
So, as we move away from the w — > oo limit where the 
spin-gapped metallic region appears for U < X, the re- 
gion becomes bounded from both below (opening of Ac) 
and above (CDW-SDW crossing point). The decrease in 
U — X at the CDW-SDW crossing may be related to the 
effective repulsion enhanced from [/ — A as studied in 
with the dynamical mean field theory (DMFT). 

If we assume, as in |0, that the Tomonaga-Luttinger 
(TL) theorvfisj apphes to the present system, in a spin- 
gapped metal the CDW correlation should behave as 
r~^p while the pairing correlation r~^/^p, so that Kp > 
1 corresponds to dominant SC correlation. While the be- 
havior of the CDW correlation with a power greater than 
unity is consistent with the results in [l^ . we also no- 




FIG. 2; Calculated powers (r;) for the correlation functions 
(oc r"") against [/ - A for (a) (A, ui) = (0.4, 5.0), (b) (A, lu) = 
(3.6,5.0) and (c) {X,uj) = (3.6,2.5) in an L = 40, half-filled 
Holstein-Hubbard chain, where U is varied along the dashed 
arrows in inset of Fig.l. Gray shadow indicates U < 0. (d) 
The result for the Hubbard model (with A = 0), where the 
curves for charge and on-site SC are degenerate. 

ticc that the power for the on-site SC is larger than that 
for CDW, which would not occur if the TL description 
persisted. Since this is observed over the entire range 
of t/ — A, this should not be due to a finite-size effect 
from the above argument. So, while a deviation from the 
TL does exist (which is not surprising since the coupling 
to phonons should make the el-el interaction effectively 
energy-dependent), the on-site SC does not become the 
most dominant correlation in the region studied here. We 
note that this robustness of CDW against SC at half- 
filled Holstein-Hubbard model has also been observed in 
a DMFT calculation for the transition temperature Q . 

About the curious near-degeneracy between the CDW 
and on-site SC correlations, we can say the following. If 
the system were purely electronic, then an electron-hole 
transformation (for down spins) maps, when the lattice is 
bipartite, the CDW correlation onto the SDW correlation 
{{S^S^}) and the on-site SC correlation onto the SDW 
correlation ((5~5+), which should be identical with 
(5^5^) due to the SU(2) symmetry). When the Hamilto- 
nian is the (half-filled) Hubbard model, the electron-hole 
transformation is the well-known attraction-repulsion 
{U ^ —U) transformation. [igj and the CDW and on- 
site pairing correlations are genuinely degenerate (as dis- 
played in FigI2tc)). Now, the system at hand is a cou- 
pled electron-phonon system, which can be mapped to 
the Hubbard model only in the limit of w — > oo. So 
the present result amounts that the degeneracy curiously 
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persists approximately for the Holstein-Hubbard model 
for finite values of u. 

Destruction of the electron-hole symmetry — This 
leads us to the following idea: a destruction of the 
electron-hole symmetry in the (one-electron) band struc- 
ture with e.g. an introduction of the second-neighbor 
hopping t' may possibly be a way to make the pairing 
correlation the single, most dominant correlation, as in 
the on-site pair correlation in the half-filled attractive 
t-t'-U model. So we introduce the second-neighbor hop- 
ping term, —t' ^(c|_,_2 ^Ci^a -\- H.c.) (= trestle lattice; 
see inset of Figl^J, in the Holstein-Hubbard model. 
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FIG. 3; (a) Correlation functions when the electron-hole 
symmetry is destroyed with t' = —0.5 for (U, A, uj) — 
(2.0,3.6,5.0). (b) Exponents against —t' for various correla- 
tion functions for (t/, A,a;) = (2.0,3.6,5.0), and (c) a similar 
plot for the attractive Hubbard model with U = —1.6. The 
inset in (b) depicts the trestle lattice. 

The result in Fig|3{a), for the same {U,X,llj) = 
(2.0,3.6,5.0) as in Fig. |2b), shows that the introduc- 
tion of t' makes the CDW correlation suppressed (i.e., 
the power increases), while the power for the on-site SC 
correlation does not significantly change. This persists, 
as shown in Fig|3Jb) plotting various exponents against 
t', up to t' — —0.5, at which the number of Fermi points 
in the non-interacting band increases from two. So in 
the region 0.25 < —t' < 0.5 the on-site SC becomes the 
single, most dominant correlation. We also note that the 
nearest-neighbor singlet pair correlation, too, becomes 
dominant (i.e., decays more slowly than the charge, spin 
or nearest-neighbor triplet pair correlations), and be- 
comes nearly as dominant as the on-site pair when —t' 
approaches 0.5. Here again, the calculated exponents 
for the Holstein-Hubbard model violates the TL relation. 



rjcDWilSC = 1, for a spin-gapped metal. 

Discussion — We have obtained an instance where su- 
perconductivity appears in a manner dependent on the 
lattice structure, while one might think that the elec- 
tronic band structure should be irrelevant for the on-site 
(an s-wave) pairing arising from the on-site U and on- 
site phonon. The result, if extensible to general dimen- 
sions, is suggestive for superconducting alkali-fuUerides 
AsCgo having an fee (i.e., electron-hole asymmetric) ar- 
ray of fuUerenes, while DMFT studies for superconduc- 
tivity primarily assume a symmetric (semielliptic) elec- 
tron density of states. Non- half- filled bands are also in- 
teresting, and the study is under way. 

This work is in part supported by a Grant-in- Aid for 
Science Research on Priority Area from the Japanese 
Ministry of Education. 



* On leave of absence from Department of Physics, Uni- 
versity of Tokyo, Hongo, Tokyo 113-0033, Japan 

[1] Z. B. Huang, W. Hanke, E. Arrigoni, and D. J. Scalapino, 
Phys. Rev. B 68, 220507(R) (2003). 

[2] M. Capone, M. Fabrizio, C. Castellani, and E. Tosatti, 
Science 296, 2364 (2002). 

[3] J. E. Han, E. Koch, and O. Gunnarsson, Phys. Rev. Lett. 
84, 1276 (2000). 

[4] J. K. Freericks and M. Jarrell, Phys. Rev. Lett. 75, 2570 
(1995). 

[5] G. S. Jeon, T.-H. Park, J. H. Han, H. C. Lee, and H.- 

Y. Choi, Phys. Rev. B 70, 125114 (2004). 
[6] O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 
[7] Y. Takada, J. Phys. Soc. Jpn. 65, 1544 (1996). 
[8] T. Hotta and Y. Takada, Physica B 230-232, 1037 

(1997) . 

[9] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[10] Y. Takada and A. Chatterjee, Phys. Rev. B 67, 

081102(R) (2003). 
[11] G. Hager, E. Jeckelmann, H. Fehske, and G. Wellein, J. 

Comp. Phys. 194, 795 (2004). 
[12] H. Fehske, G. Wellein, G. Hager, A. Weifie, and A. R. 

Bishop, Phys. Rev. B 69, 165115 (2004). 
[13] E. Jeckelmann and S. R. White, Phys. Rev. B 57, 6376 

(1998) . 

[14] M. Tezuka, R. Arita, and H. Aoki, Physica B 359-361, 
708 (2005). 

[15] J. Solyom, Adv. Phys. 28, 201 (1979); V. J. Emery 
in Highly Conducting One-Dimensional Solids, ed. by 
J. T. Devreese et al. (Plenum, New York, 1979), p. 247; 
H. J. Schulz in Mesoscopic Quantum Physics, ed. by 
E. Akkermans et al. (Elsevier, Amsterdam, 1995), p. 533. 

[16] K. Kuroki, T. Kimura, and H. Aoki, Phys. Rev. B 54, 
R15641 (1996). 

[17] G. Sangiovanni, M. Capone, C. Castellani, and M. Grilli, 

Phys. Rev. Lett. 94, 026401 (2005). 
[18] R. T. Clay and R. P. Hardikar, Phys. Rev. Lett. 95, 

096401 (2005). 
[19] Y. Nagaoka, Prog. Theor. Phys. 52, 1716 (1974). 



